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Abstract: We study bulk fermion fields in the localized gravity model with 
non-factorizable metric recently proposed by Randall and Sundrum, and Gog- 
berashvili. In addition to a tower of weak-scale Kaluza-Klein states we find a zero 
mode for any value of the fundamental fermion mass. If the fermion mass is larger 
than half the curvature of the compact dimension, the zero mode can be localized 
on the "hidden" 3-brane in the Randall-Sundrum model. Identifying this mode 
with a right-handed neutrino provides a new way for obtaining small Dirac neutrino 
masses without invoking a see-saw mechanism. Cancellation of the parity anomaly 
requires introducing an even number of bulk fermions. This naturally leads to a 
strong hierarchy of neutrino masses and generically large mixing angles. 
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1. Introduction 

Theories with extra spatial dimensions have received great attention recently, when 
it was shown that they could provide a solution to the gauge-hierarchy problem. 
If space-time is a product of Minkowski space with n compact dimensions, with 
Standard Model fields localized in the three extended spatial dimensions (i.e., on a 
3-brane) and gravity propagating in the extra space, then the strength of gravity 
on the 3-brane is governed by an effective Planck scale Mp { = M n+2 V n , where M 
is the fundamental scale of gravity and V n the volume of the compact space [1]. If 
this space is sufficiently large, the fundamental scale M can be of order 1 TeV, thus 
removing the large disparity between the gravitational and the electroweak scales. 

An intriguing alternative to the above scenario invokes a non-factorizable geom- 
etry with a metric that depends on the coordinates of the extra dimensions [2, 3]. In 
the simplest scenario due to Randall and Sundrum (RS) one considers a single extra 
dimension, taken to be a S 1 /Z 2 orbifold parameterized by a coordinate y = r c 0, 
with r c the radius of the compact dimension, — n < <fi < ir, and the points (x, 0) and 
(x,—(j)) identified [3]. There are two 3-branes located at the orbifold fixed points: 
a "visible" brane at <fi = it containing the Standard Model fields, and a "hidden" 
brane at = 0. The solution of Einstein's equations for this geometry leads to the 
non-factorizable metric 

ds 2 = e - 2fo ^l %u dx^dx" - r\ d0 2 , (1.1) 

where x M are the coordinates on the four-dimensional surfaces of constant 0, and the 
parameter k is of order the fundamental Planck scale M. (This solution can only 
be trusted if k < M, so the bulk curvature is small compared with the fundamental 
Planck scale.) The two 3-branes carry vacuum energies tuned such that Kds = 
—Vhid = —24M 3 k, which is required to obtain a solution respecting four- dimensional 
Poincare invariance. In between the two branes is a slice of AdS 5 space. 

With this setup, the effective Planck scale seen by particles confined to four- 
dimensional space-time is Mpj = (M 3 //c)(l — e~ 2krcW ), which is of order the funda- 
mental scale M. Unlike the scenarios with large extra dimensions considered in [1], 
the scale M is therefore not of order the weak scale. However, the "warp factor" 
e -2kr c \4>\ j n me t r j c (1.1) has important implications for the masses of particles 
confined to the visible brane. The Lagrangian for these particles depends on the 
induced metric on the brane, g™ = e~ 2krc7r r/^, and after field renormalization any 
mass parameter m in the fundamental theory is promoted into an effective mass pa- 
rameter m = e~ kTcn m governing the physical properties of particles on the brane [3]. 
With kr c pa 12 this mechanism produces weak-scale physical masses and couplings 
from fundamental masses and couplings of order the Planck scale. As a consequence 
of the warp factor, the Kaluza-Klein excitations of gravitons have weak-scale mass 
splittings and couplings [4, 5]. This is in contrast with the Kaluza-Klein spectrum 
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of gravitons propagating in large extra dimensions, which consists of a large number 
of light modes (with splittings of order the compactification scale) with gravitational 
couplings. The same properties (i.e., weak-scale masses and couplings) are shared 
by bulk scalars and vector particles propagating in the extra dimension [6, 7, 8]. 

In order for this model to provide a viable solution to the hierarchy problem it 
is important to address the question of how to stabilize the radius r c of the extra 
dimension, and the related question of the potentially disastrous cosmology of a 
visible universe confined to a brane with negative tension [9, 10, 11] (see also [12, 
13]). A mechanism for radius stabilization utilizing a bulk scalar field has been 
proposed in [14]. In the presence of such a mechanism, one finds a conventional 
cosmological expansion for temperatures below the weak scale [15]. The couplings 
of the radion field to Standard Model particles may have interesting implications for 
collider searches [15, 16]. Other phenomenological consequences of the RS model 
have been discussed in [17], and an alternative model avoiding the negative-tension 
brane has been proposed in [5]. 

Resolving the hierarchy problem by introducing extra dimensions poses new 
challenges. For instance, operators mediating proton decay, lepton-number violation 
or flavor- changing neutral currents must be sufficiently suppressed. Likewise, the see- 
saw mechanism for generating small neutrino masses cannot be invoked if the highest 
energy scale governing physics on the visible brane is the weak scale. However, there 
is now increasing evidence that the atmospheric neutrino anomaly [18, 19, 20] and 
the solar neutrino problem [21, 22] are explained in terms of neutrino oscillations, 
which require small but non-vanishing neutrino masses. Several four- dimensional 
alternatives to the see-saw mechanism not requiring a high-energy scale have been 
proposed, such as radiatively generated neutrino masses [23] and composite models 
[24]. However, it would be interesting to find new mechanisms that are intrinsically 
higher dimensional. In the context of models with large extra dimensions ideas in this 
direction have been presented in [25, 26], and some concrete models have been worked 
out in [27, 28, 29]. They contain a massless Standard Model singlet propagating in the 
bulk of the extra compact space, which serves as a right-handed neutrino. Then the 
effective four-dimensional Yukawa coupling is suppressed by a volume factor 1 / \/Vn, 
reflecting the small overlap between the right-handed neutrino in the bulk and the 
left-handed one on the 3-brane. By construction, this factor provides a suppression 
of neutrino masses of order v/Mp\, reminiscent of the see-saw mechanism. However, 
this idea does not work in a scenario with small extra dimensions such as the RS 
model, simply because of the lack of a volume suppression factor. 

In this letter we investigate the possibility of incorporating bulk fermions in 
the RS model. As in the case of scalars [6] and vector fields [7, 8] propagating 
in the compact dimension we find that the Kaluza-Klein modes have weak-scale 
masses even though the fermion mass in the fundamental, five-dimensional theory 
is of order the Planck scale. The fermion case is more interesting, however, because 
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the extension of the Dirac algebra to five dimensions leads to a different propagation 
of left- and right-handed modes. After imposing the orbifold boundary conditions 
the geometry supports a left-handed or a right-handed zero mode for any value of 
the fundamental fermion mass, one of which can be localized on the hidden brane of 
the RS model. This is different from the scalar and vector cases, where zero modes 
exist only for vanishing mass in the fundamental theory. The localization of a right- 
handed zero mode on the hidden brane provides a new mechanism for obtaining small 
neutrino masses, which can be realized by coupling the Higgs and left-handed lepton 
fields of the Standard Model, localized on the visible brane, to a right-handed fermion 
in the bulk. The neutrino mass can be tuned over many orders of magnitude by a 
small change of the bulk fermion mass. Moreover, cancellation of the parity anomaly 
[30, 31] forces us to introduce an even number of bulk fermions. This naturally leads 
to a neutrino mass hierarchy and potentially large mixing angles. 



2. Bulk fermions 

Our starting point is the action for a Dirac fermion with mass m of order the funda- 
mental scale M propagating in a five-dimensional space with the metric (1.1), which 
we write in the form 



J d 4 x J d<f) 



G\E t 



A 



— m sgn(0) 



(2.1) 

where G = det (Gab) = r^e~ 8cr with a = kr c \<f)\ is the determinant of the metric. 
We use capital indices A,B, . . . for objects defined in curved space, and lower-case 
indices a, b, . . . for objects defined in the tangent frame. The matrices 7 a = (7^, 175) 
provide a four-dimensional representation of the Dirac matrices in five- dimensional 
flat space. The quantity E£ = diag(e CT , e CT , e CT , e CT , l/r c ) is the inverse vielbein, and 
UbcA is the spin connection. Because in our case the metric is diagonal, the only non- 
vanishing entries of the spin connection have b = A or c = A, giving no contribution 
to the action in (2.1). 

The sign change of the mass term under — > —0 is necessary in order to conserve 
0-parity, as required by the Z2 orbifold symmetry of the RS model. Such a mass term 
can be obtained, e.g., by coupling the fermion to a pseudoscalar (under 0-parity) bulk 
Higgs field. For a single bulk fermion in five dimensions 0-parity is broken at the 
quantum level, giving rise to the so-called parity anomaly [30, 31]. To cancel this 
anomaly, we will later consider an even number of fermion fields. 



x We do not include a five-dimensional Majorana mass term of the form ^> T C ^> in the action, 
because later we will assign lepton number to the bulk fermion. 
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Using an integration by parts, and defining left- and right-handed spinors ^l,r = 
|(1 =f 75)^/, the action can be written as 

S = Jd 4 xJ d0r c je- 3CT (* L i0tt L + - e- 4CT msgn(0) (tt L tt fl + 

"^l (e-^dj, + e- 4CT ) - V R (V 4 ^ + e" 4CT ) tt L ] } , (2.2) 



1 

2^ 



where we impose periodic boundary conditions ^l,r{x,h) = ^l,r(x, — n) on the 
fields. The action is even under the Z 2 orbifold symmetry if ^ L (x,(f)) is an odd 
function of <p and ^r(x, 0) is even, or vice versa. To perform the Kaluza-Klein 
decomposition we write 



e 2a 



<t>) = E ti' R (x) ~~i=- fn' R (<P) • (2-3) 



Because of the Z2 symmetry of the action it is sufficient to restrict the integration 
over (f) from to ir. The behavior of the solutions for negative is then determined 
by their Z 2 parity. {/^ (</>)} and {/,f (</>)} are two complete, orthonormal (with a 
scalar product defined below) sets of functions on the interval G [0,7r], subject 
to certain boundary conditions. We will construct them as the eigenfunctions of 
hermitian operators on this interval. Inserting the ansatz (2.3) into the action and 
requiring that the result take the form of the usual Dirac action for massive fermions 
in four dimensions, 

S = Y1 d 4 x {^ n (x) ip ij n (x) -m n 4> n (x)4> n (x)} , (2.4) 

where ip = ipi + 'ipR (except for possible chiral modes) and m n > , we find that the 
functions fn' R {4>) must obey the conditions 

J d0eV^*(0) /„ L (0) = / d0e CT /^(0) = 5 mn , 



(±i d+ - m) ft R {4>) = ~rn n e<7,?' L (0) . (2.5) 

The boundary conditions /^*(0) /,f (0) = f^i^) fni 71 ) — 0> which follow since either 
all left-handed or all right-handed functions are Z 2 -odd, ensure that the differential 
operators (±^- d^ — m) are hermitian and their eigenvalues m n real. (Since the equa- 
tions are real, the functions fn ,R {4>) could be chosen real without loss of generality.) 

It is convenient to introduce the new variable t = ee a G [e, 1] with e = e~ kTc7T , 
rescale fn' R ((f>) — > Vkr c e f^' R (t), and define the quantities 

m m n m n k 

U = I> Xn = 7k=T e ■ (2 ' 6) 
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The small parameter e ~ 10 16 sets the ratio between the electroweak and the 
gravitational scales. The two conditions in (2.5) now become 

i i 

/ & C(t) ft(t) = jdt C(t) /* (t) = 6 mn , 

e e 

(±td t - v)ft\t) = -x n tf R > L (t) , (2.7) 

and the boundary conditions are fm*( e ) fn( e ) = /m*(l) /?f(l) = 0- The system of 
coupled, first-order differential equations for fn' R (t) implies the second-order equa- 
tions 

~t 2 d 2 t + xlt 2 - v{v =F 1)] ft R (t) = . (2.8) 

Dimensional analysis shows that the eigenvalues x n are of order unity, corresponding 
to weak-scale fermion masses m n in the four- dimensional theory. 

The solution of the differential equations is straightforward. We start by looking 
for zero modes, i.e., solutions with x n = 0. In this case the first-order equations in 
(2.7) decouple. The properly normalized solutions are 

fo' R (t) = fo' R (l) t ±u oc e^'W , |/ ^(1)| 2 = ^±§^ . (2.9) 

Since these are even functions of 0, which do not vanish at the orbifold fixed points, 
only one of the zero modes is allowed by the orbifold symmetry. This mode exists 
irrespective of the value of the fermion mass m in the five-dimensional theory. Note 
that for v > \ the right-handed zero mode has a very small wave function on the 
visible brane: f R (l) oc e u ~?. This property will allow us to obtain small neutrino 
masses. The presence of fermion zero modes should not come as a surprise, since it is 
well known that in flat space-time they are associated with domain walls [33]. In our 
model the domain walls are provided by the 3-branes of the RS model, which separate 
the regions with a different sign of the fermion mass term. The functions Pf,/^^) 
and P R f^(t), with P LyR = |(1 =F 75), can be associated with the "fermionic" and 
"bosonic" degrees of freedom of a supersymmetric, quantum-mechanical system [34]. 
The supersymmetry generators are Q — (d t — v /i)7 Pj, and = —(d t + 17^)7° Pr, 
and the Kaluza-Klein modes are the eigenstates of the Hamiltonian {Q,Q^}- This 
explains why left- and right-handed modes have the same eigenvalues x n . The two 
zero modes correspond to the ground-state solutions of the supersymmetric Hamil- 
tonian. In our case, the requirement of orbifold symmetry allows only one of these 
solutions to be present. 

The solutions of the differential equations (2.8) for the case x n > are Bessel 
functions. For convenience we assume that v 7^ | + iV with an integer N. Then the 
most general solutions can be written in the form 

fn' R (t) = Vt [a L n > R Ji TU (x n t) + b L n < R J_ h±v {x n t)] . (2.10) 
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Option 


Eigenvalues 


X n > 


a n 




Zero Mode 


L, v< \ 




= 







R 


h,v>\ 


^v-\ ( X n) 


= 


' 




R 


R 




= 







L 



Table 1: Bulk fermion solutions for the two choices of boundary conditions, in the limit 
where e = e~ kTc7T — > 0. The wave functions fn' R (t) take the form (2.11) with coefficients 
a^ ,R given in the third and fourth columns. 



For the special values v = \ + N the solutions are superpositions of Bessel functions 
of the first and second kind, which can obtained from our results using a limiting 
procedure. The two functions fn(t) and f R (t) are not independent, since they are 
coupled by the first-order differential equations in (2.7), which imply = a R and 
b R = — a^. Hence, the solutions take the form 

fn (*) = [a% Ji_ v (x n t) + a R J_i +V (x n t)] , 

f R (t) = Vt [a R Ji 2+U (x n t) - a L n J_ k _ v (x n t)] . (2.11) 

To proceed we must specify the boundary conditions at the locations of the 3- 
branes. This will give rise to a discrete spectrum of Kaluza-Klein modes. Orbifold 
symmetry allows two choices of boundary conditions: either all left-handed fields 
are odd under 0-parity and all right-handed ones even ("option L"), or all right- 
handed fields are odd and all left-handed ones even ("option R"). In the first case 
the boundary conditions are /^(e) = /^(l) = 0, and in the second case f R (e) = 
f R (l) = 0. Which of these choices is realized in nature is a question that cannot be 
answered without understanding the physics on the 3-branes, which is beyond the 
scope of the field-theory model suggested in [3]. The two cases are straightforward to 
analyze. Using the asymptotic behavior of the Bessel functions, J n (x) ~ x n as x — > 0, 
it follows that in the limit e — > only one of the two terms in the wave functions 
in (2.11) remains. Taking e — > is an excellent approximation unless we were to 
consider integrals of the functions fn' R (t) with weight functions that are singular 
as t — > 0. Table 1 summarizes the results for the eigenvalues and eigenfunctions 
of the non-zero modes in that limit. The solutions shown correspond to positive 
and must be extended to negative in accordance with the orbifold symmetry. For 
option L the results take a different form depending on whether v < \ or v > |, 
as indicated. The second column in the table shows the equation that determines 
the eigenvalues x n . In the next two columns we give the values of the coefficients 
a^' R of the properly normalized solutions. The normalization constants Af a (x) obey 
\Af a (x)\ 2 = 2/[J a (x)} 2 . The last column shows the chirality of the zero mode. The 
zero-mode wave functions can be recovered by taking the limit x n — > 0; however, the 
normalization constants do not apply in this case. In figure 1 we show the first few 
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Figure 1: Right-handed (solid) and left-handed (dashed) Kaluza-Klein modes with n < 2 
for v = m/k = 0.45 (left) and 0.55 (right), with boundary conditions such that all left- 
handed fields vanish at t = e and t = 1 (option L). We show exact results obtained with 
e = 10" 16 . In both cases, the corresponding eigenvalues are x\ ~ 2.49 and £2 ~ 5.60. 



Kaluza-Klein modes for option L and two values of the parameter v just below or 
above the critical value v — |. The important point to notice is the localization of 
the right-handed zero mode f^(t) on the hidden brane (at t — e) for v > \. 

For the special case of integer v the exact solutions for the wave functions can 
be expressed in terms of trigonometric functions. As an example, we quote results 
for v = and v = 1 choosing for the boundary conditions option L, which will be 
of special relevance to our discussion below. In both cases the non-zero eigenvalues 
are given by x n = mr/(l — e) with an integer n > 1, and the left-handed solutions 
are f^(t) = J\fsm[x n (t — e)], where \Af\ 2 = 2/(1 — e). For v — 0, the right-handed 
solutions are given by 

fo R (t) = ^ , f*(t) = -Afcos[x n (t - e)] , (2.12) 

whereas for v — 1 they take the form 

/*(*) = ^ , /«(r) = AT ( 5i "'^- g " - C08 M< - e) l) ■ (2.13) 

3. Yukawa interactions and neutrino phenomenology 

We will now show how including a sterile bulk fermion in the RS model can provide 
a mechanism for obtaining small Dirac neutrino masses, which is quite different from 
the see-saw mechanism. We focus first on a single fermion generation and consider a 
scenario where all matter and gauge fields charged under the Standard Model gauge 
group are confined to the visible brane at <fi = tt, whereas a gauge-singlet fermion 
field propagates in the bulk. After integration over the compact extra dimension we 
obtain a tower of four-dimensional Kaluza-Klein fermions in the four-dimensional 
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theory, as shown in (2.4). We choose boundary conditions such that there is a right- 
handed zero mode (option L) with wave function f^(t) given in (2.9). Only this 
choice will lead to an interesting neutrino phenomenology. 

Omitting gauge interactions, the action for a Higgs doublet H = (0i, ^2), a left- 
handed lepton doublet L = {y L ,e L ) and a right-handed lepton e R localized on the 
visible brane is 

S = |d 4 x v ^-{c^o t ^o-A(|i/o| 2 -^) 2 } 

+ / d 4 x V-fiVis {ZoT^Lo + zmYd^m - (y e L H e R0 + h.c.) } , (3.1) 

where g^ s = e 2kTc7T rf v is the induced metric on the brane, ^—gv\s = det(— g™) = 
e ~4kr c n^ an j _ fiv^ — 7r) 7° = e kTc7T 7 M . To restore a canonical normalization of 
the fields on the brane we must perform the rescalings H — > e kr ' c7T H, L — > ei kTcn L 
and — > e^ kr ' c7T e R , upon which the action takes the form 

S = jd 4 x [d^d'H - X (\H\ 2 - v 2 ) 2 } 

+ J d 4 x {LiftL + e R ip e R - (y e LHe R + h.c.) } , (3.2) 

where v = e~ kr ' c7T v . The remarkable feature noted in [3] is that all dimensionful 
parameters such as the Higgs vacuum expectation value get rescaled by the warp 
factor and turned from Planck-scale into weak-scale couplings, whereas dimensionless 
parameters such as A and y e remain unchanged. 

We now introduce a Yukawa coupling of the bulk fermion with the Higgs and 
lepton fields. With our choice of boundary conditions all left-handed Kaluza-Klein 
modes vanish at the visible brane, so only the right-handed modes can couple to the 
Standard Model fields on the brane. However, in a more realistic scenario which takes 
into account a finite width of the 3-branes there will most likely be a non-zero (and 
indeed sizeable) overlap of the left-handed modes with the Standard Model fields. 
Hence, in order to avoid weak-scale neutrino masses and lepton-number violating 
interactions we assign lepton number L — 1 to the bulk fermion state. Then the only 
gauge-invariant coupling is of the form 

S Y = - Jd 4 xV^ s {Y 5 L (x)H (x)^ R (x,TT)+h.c.} , (3.3) 

where H = ia 2 H*, and the Yukawa coupling Y 5 is naturally of order M~ 1//2 , with M 
the fundamental Planck scale of the theory. Rescaling the Standard Model fields in 
the way described above, and inserting for the bulk fermion the Kaluza-Klein ansatz 
(2.3), we find 

S y = - E [ d ' x {y n L(x)H(x)^(x) + h.c.} (3.4) 

n>0 J 
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with the effective Yukawa couplings 

yn = VkY s tf(l) = Y B f?(l), (3.5) 

where Y 5 is naturally of order unity. After electroweak symmetry breaking, this 
Yukawa interaction gives rise to a neutrino mass term Ai ip R + h.c, which in the 
basis ip u L = [vli ipi, ■ ■ ■ , VnO an< ^ = {^oi ^11 ■ ■ ■ i V'n )■> w ith n — > oo, takes the form 



mi ... 

: '•• 
V ... m n J 



(3.6) 



As a consequence, there will be a mixing of the Standard Model neutrino ul with the 
heavy, sterile (with respect to the Standard Model gauge interactions) bulk neutrinos 
■ipn- The Kaluza-Klein excitations of the bulk fermion have masses m n of order the 
weak scale v. Thus, in order to obtain a light neutrino we need | j/ 1 ^ 1) which 
requires having a very small wave function of the zero mode on the visible brane, 
i.e., |/(f(l)| <C 1. But this is precisely what happens if the fundamental fermion 
mass m satisfies the condition m > k/2. Since, as mentioned earlier, the curvature k 
must be smaller than the fundamental scale M, this is a natural requirement in the 
context of the RS model. 

In order to study the properties of the physical neutrino states we diagonalize 
the squared mass matrix The eigenvalues of this matrix are the squares of 

the physical neutrino masses, and the unitary matrix U defined such that WAiAi^U 
is diagonal determines the left-handed neutrino mass eigenstates via i\) v L = C/^£ hys . 
We denote by m v the mass of the lightest neutrino z/£ hys and define a mixing angle 
9 U such that v L = cos 6^ i/£ hys + . . ., where the dots represent the admixture of heavy, 
sterile bulk states. To leading order in the small parameter |y | C 1 we obtain 

m u — v\yo\ cos9 u , tan 2 6 v — "Y] . (3.7) 

Since the mixing angle is constrained by experiment to be very small (see below), it 
follows from (2.9) and (3.5) that 

m u ~ V2u - 1 \Y 5 \ e u ~^ v ~ M (^-Y +2 ; u>-. (3.8) 

This result is remarkable, as it provides a parametric dependence of the neutrino 
mass on the ratio of the electroweak and Planck scales that is different from the 
see-saw relation m v ~ v 2 /M, except for the special case where v — |. This flexibility 
allows us to reproduce a wide range of neutrino masses without any fine tuning. 
For instance, taking v/M = 10~ 16 , the phenomenologically interesting range of m u 
between 10 _5 eV and 10 eV can be covered by varying v between 1.1 and 1.5. 
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The measurement of the invisible width of the Z Q boson, which yields n u = 
2.985 ±0.008 for the apparent number of light neutrinos [35], implies that the mixing 
angle 6 V must be of order a few percent. For instance, assuming an equal admixture of 
sterile neutrinos for the three generations of light neutrinos, we obtain n v = 3 cos 2 u 
and hence tan 2 6 u = 0.005 ± 0.003. From table 1 it follows that with our choice 
of boundary conditions the wave functions of all excited right-handed Kaluza-Klein 
modes obey (1)| = V2 (for e — > 0). We thus obtain 




actly and yields l/(2v + 1). To satisfy the bound on the mixing angle for v = 0(1) 
requires that v \Y 5 \/k < 0.1, which is possible without much fine tuning. We em- 
phasize, however, that it would be unnatural to have the dimensionless combination 
volY^l/k much less than unity, so a mixing angle 6 U not much smaller than the cur- 
rent experimental bound is a generic feature of our scenario, which can be tested by 
future precision measurements. 

So far we have shown how a right-handed bulk fermion can give a small Dirac 
mass to a Standard Model neutrino. We now generalize this mechanism to three 
neutrino flavors and more than one bulk fermion. Interestingly, such a generalization 
is forced upon us by the requirement that the parity anomaly for fermions in an 
odd number of dimensions vanish. When an odd number of bulk fermions in five 
dimensions are coupled to a gauge field or gravity, the 0-parity of the action (2.1) is 
broken at the quantum level [30, 31]. To obtain a minimal model that is anomaly free 
we thus introduce two bulk fermions with boundary conditions option L, so there 
are two massless right-handed zero modes. 2 In order to explain the atmospheric 
and solar neutrino anomalies in terms of neutrino oscillations one needs two very 
different mass-squared differences: Am^ <C Am^, where Amy = to 2 ,. — m 2 _., and by 
convention m vi < m V2 < m U3 . This requires a minimum of two massive neutrinos; 
however, the third neutrino can be massless. In our minimal model this is indeed 
what happens. Although it is perhaps unconventional to consider a scenario where 
the number of right-handed neutrinos does not match the number of left-handed ones, 
we will see that our model explains successfully the known features of the neutrino 
mass and mixing parameters. 

In order to explore this minimal model in more detail we ignore, for simplicity, 
the heavy Kaluza-Klein excitations of the bulk fermions and focus only on the zero 
modes. As mentioned above, the admixture of weak-scale sterile neutrino states 

2 More complicated models with four or more bulk states are possible. These states could be 
subject to different boundary conditions. If we impose lepton number, only the right-handed modes 
can couple to the Standard Model fields. At least two right-handed zero modes are needed for a 
successful neutrino phenomenology. 
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must be strongly suppressed. It is natural to allow for the possibility that the two 
bulk fermions have different masses m\ > m 2 (of order the Planck scale) in the 
fundamental theory, and that they couple with similar strength to the three left- 
handed neutrino flavors. According to (2.9) and (3.5), the effective Yukawa couplings 
of the two right-handed zero modes ifj^' 1 and ?/>(f' 2 can then be parameterized as 
Xf€ Ul ~^ and yje u ' 2 ~^ with z/j = rrii/k (for i — 1,2) and flavor-dependent couplings 
Xf, i/f (with / = e,/i, r) of order unity. Note that the Yukawa couplings of the 
two zero modes have a very different magnitude: Xf/yf = 0(e Ul ~ u ' 2 ). The resulting 
neutrino mass term Tp v L M.ip u R + h.c. in the truncated basis ip v L = (z/f, v^) and 

n (<■■!!". c!;- 2 ) is 

Ve \ 



M = ve U2 ~* | e^x^ Vil 

VrJ 



(3.10) 



Diagonalizing the matrix M.M.^ to leading order in e we find that the physical neu- 
trino mass eigenstates comprise a massless left-handed neutrino vi, a very light Dirac 
neutrino with mass squared 

\uX- + W + \vA 2 \m) ' (J11 > 

and a light Dirac neutrino with mass squared 

^e 2 -- 1 (|y e | 2 + M 2 + \y T \ 2 ) ~ M 2 (^£f • (3.12) 



In (3.11) we use the short-hand notation [ij] = x^y^ — Xjyi. Since the lightest 
neutrino is massless it follows that Am^ = m 2 2 and Am| 2 ~ m 2 ^, and the ratio 
Ama/Am^ ~ (t> /M) 2 ^ 1 "^. An interpretation of the solar neutrino anomaly in 
terms of neutrino oscillations based on the MSW effect [36] yields values of Am^ 
in the range 10~ 6 -10~ 5 eV 2 , whereas oscillations in vacuum would require a smaller 
value of order 10 _10 eV 2 [22]. Such masses can be reproduced in our model by set- 
ting Vi m 1.34-1.37 and v\ ~ 1.5, respectively. An explanation of the atmospheric 
neutrino anomaly in terms of neutrino oscillations yields Am 2 2 in the range 5 • 10 -4 - 
6 ■ 10~ 3 eV 2 [21], which we can reproduce by taking v 2 ~ 1.27-1.29. In other words, 
we can understand the observed hierarchy of the experimentally favored neutrino 
masses in terms of a small difference of the bulk fermion masses in the fundamental 
theory. Note that in our minimal model the lightest neutrino is massless. This can 
be changed by introducing four (or more) bulk fermion states with more than two 
right-handed zero modes, in which case also the lightest neutrino becomes massive, 
with m 2 Vl <ti m 2 2 . 

Despite the fact that a strong neutrino mass hierarchy is a generic feature of 
our model, the mixing matrix U relating the neutrino flavor and mass eigenstates 
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does not contain any small parameter. Defining v$ = Y^=i Ufi v-i we find that all the 
entries U fi are of order unity. In the limit e^Owe obtain 



/u el u e2 u e3 \ 



u 



u, 



IJ.2 



u, 



,i3 



\u Tl u t2 U t3 J 





[e(i]-y* M 




JVi7V 2 


[re]* 


v* [^ T l -vt M 


Ni 


NiN 2 


M* 


yi[re]-y*\Mr] 


Ni 


NiN 2 



N 2 \ 

Mr. 

N 2 

Vr I 

N 2 / 



(3.13) 



where Nf = |[e/i]| 2 + | [/xt] | 2 + | [-re] | 2 and iVf = |y e | 2 + |^| 2 + \y T \ 2 - A mixing matrix of 
this type, which lacks the strong hierarchy of the quark mixing matrix, can account 
for the experimental constraints on the neutrino mixing angles. (In fact, it has been 
pointed out that a fair fraction of random Dirac mixing matrices is consistent with 
these constraints [37].) The precise form of these constraints depends on how the data 
are analyzed, and whether the solar and atmospheric neutrino anomalies individually 
are interpreted in terms of two-neutrino or three-neutrino mixing. Constraints from 
the CHOOZ reactor experiment [38] combined with the atmospheric neutrino data 
imply that \U e3 \ 2 ^ few % [39], which means that \y e \ should be less than and 



In the limit where \y e \ 2 <C |y M | 2 + \y T \ 2 , the mixing angles &12 and #23 responsible 



for v e <-> and <-> v T oscillations obey the approximate relations 



sin 2 2#i2 ~ 
sin 2 2#23 ~ 



4| 


X e 


2 (M 2 + I 


Vr 


2 )IMI 2 


[\x e 


\ 2 iW + 


Ur 


l 2 ) + IMI 2 ] 2 



4|^| 2 |l/r| 



2^2 • 



(3.14) 



The atmospheric neutrino anomaly is best explained by near-maximal <-> v T mix- 
ing, such that sin 2 26 l 2 3 > 0.82 [20]. This implies 0.64 < \y^/y T \ < 1-57, which 
clearly does not pose a problem for our model. Likewise, a large-mixing-angle 
solution to the solar neutrino problem requires sin 2 2#i2 > 0.75 [22], which yields 
0.58 < l^el^l^l 2 + |y T | 2 /| [yur] I < 1.73. The small-mixing-angle MSW solution, on 
the other hand, prefers sin 2 20 12 
2 and \x^y T 



\Xe\\J\Vn\ + ml 
achieved by having \x e 



10~ 2 , which would require that the quantities 
x T y^\ differ by about a factor 20. This could either be 
\x^ T \, or via a near degeneracy of x^y T and Xry^. 



4. Conclusions 

We have studied bulk fermion solutions in the localized gravity model with non- 
factorizable geometry introduced by Randall and Sundrum to solve the gauge-hierar- 
chy problem. Similar to the case of scalar or vector fields propagating in the extra 
compact dimension, we have found that the Kaluza-Klein modes have weak-scale 
masses even if the fermion mass in the fundamental, five- dimensional theory is of 
order the Planck scale. However, a distinct feature of bulk fermion solutions is the 
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possible presence of zero modes due to the fact that the 3-branes in the Randall- 
Sundrum model act as domain walls. 

Our most important finding is that, if the fundamental mass m in the five- 
dimensional theory is larger than half the curvature k of the compact space, an 
appropriate choice of the orbifold boundary conditions leads to a right-handed zero 
mode localized on the hidden brane, whose wave function on the visible brane is 
power-suppressed in the ratio of the weak scale to the fundamental Planck scale. 
Coupling the Higgs and left-handed lepton fields of the Standard Model, localized 
on the visible brane, with a bulk right-handed neutrino provides a new mechanism 
for obtaining small neutrino masses. Remarkably, this mechanism leads to a gener- 
alization of the see-saw formula with a different parametric dependence on the ratio 
v/M, which can easily reproduce neutrino masses in the range 10 _5 eV to 10 eV. 
Without much fine tuning the mixing of the Standard Model left-handed neutrino 
with sterile, weak-scale Kaluza-Klein excitations of the bulk fermion can be made 
consistent with experimental bounds. However, a generic prediction of our model is 
that such a mixing should occur at a level not much below the present bound. 

Finally, we have shown that with an even number of bulk fermions one can obtain 
viable models of neutrino flavor oscillations, which naturally predict a mass hierarchy 
and a neutrino mixing matrix not containing any small parameter. A minimal imple- 
mentation of this scenario consists of two right-handed neutrinos, identified with the 
zero modes of two bulk fermions with slightly different masses in the five- dimensional 
theory, coupled to the three left-handed neutrinos of the Standard Model. In this 
model we obtain a massless left-handed neutrino and two massive Dirac neutrinos 
with a large mass hierarchy and generically large mixing angles. 
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